Large scale phenotyping of long COVID inflammation reveals mechanistic subtypes of disease (preprint)

One in ten SARS-CoV-2 infections result in prolonged symptoms termed "long COVID", yet disease phenotypes and mechanisms are poorly understood. We studied the blood proteome of 719 adults, grouped by long COVID symptoms. Elevated markers of monocytic inflammation and complement activation were associated with increased likelihood of all symptoms. Elevated IL1R2, MATN2 and COLEC12 associated with cardiorespiratory symptoms, fatigue, and anxiety/depression, while elevated MATN2 and DPP10 associated with gastrointestinal (GI) symptoms, and elevated C1QA was associated with cognitive impairment (the proteome of those with cognitive impairment and GI symptoms being most distinct). Markers of neuroinflammation distinguished cognitive impairment whilst elevated SCG3, indicative of brain-gut axis disturbance, distinguished those with GI symptoms. Women had a higher incidence of long COVID and higher inflammatory markers. Symptoms did not associate with respiratory inflammation or persistent virus in sputum. Thus, persistent inflammation is evident in long COVID, distinct profiles being associated with specific symptoms.

Empagliflozin in patients admitted to hospital with COVID-19 (RECOVERY): a randomised, controlled, open-label, platform trial (preprint)

Background: Empagliflozin has been proposed as a treatment for COVID-19 on the basis of its anti-inflammatory, metabolic and haemodynamic effects. Methods: In this randomised, controlled, open-label trial, several possible treatments are compared with usual care in patients hospitalised with COVID-19. Eligible and consenting adults were randomly allocated in a 1:1 ratio to either usual standard of care alone or usual standard of care plus empagliflozin 10mg once daily for 28 days or until discharge using web-based simple (unstratified) randomisation with allocation concealment. The primary outcome was 28-day mortality. On 3 March the independent data monitoring committee recommended that the investigators review the data and recruitment was consequently stopped on 7 March. The trial is registered with ISRCTN (50189673) and clinicaltrials.gov (NCT04381936). Findings: Between 8 July 2021 and 6 March 2023, 4271 patients were randomly allocated to receive either empagliflozin (2113 patients) or usual care alone (2158 patients). Overall, 289 (14%) patients allocated to empagliflozin and 307 (14%) patients allocated to usual care died within 28 days (rate ratio 0.96; 95% confidence interval [CI] 0.82-1.13; p=0.64). There was no evidence of significant differences in duration of hospitalisation (median 8 days vs. 8 days) or the proportion of patients discharged from hospital alive within 28 days (79% vs. 78%; rate ratio 1.03; 95% CI 0.96-1.10; p=0.44). Among those not on invasive mechanical ventilation at baseline, there was no evidence of a significant difference in the proportion meeting the composite endpoint of invasive mechanical ventilation or death (16% vs. 17%; risk ratio 0.95; 95% CI 0.84-1.08; p=0.44). Interpretation: In adults hospitalised with COVID-19, empagliflozin was not associated with reductions in 28-day mortality, duration of hospital stay, or risk of progressing to invasive mechanical ventilation or death.

Higher dose corticosteroids in hospitalised COVID-19 patients with hypoxia but not requiring ventilatory support (RECOVERY): a randomised, controlled, open-label, platform trial (preprint)

Background: Low-dose corticosteroids have been shown to reduce mortality for hypoxic COVID-19 patients requiring oxygen or ventilatory support (non-invasive mechanical ventilation, invasive mechanical ventilation or extra-corporeal membrane oxygenation). We evaluated the use of a higher dose of corticosteroids in this patient group. Methods: This randomised, controlled, open-label platform trial (Randomised Evaluation of COVID-19 Therapy [RECOVERY]) is assessing multiple possible treatments in patients hospitalised for COVID-19. Eligible and consenting adult patients with clinical evidence of hypoxia (i.e. receiving oxygen or with oxygen saturation <92% on room air) were randomly allocated (1:1) to either usual care with higher dose corticosteroids (dexamethasone 20 mg once daily for 5 days followed by 10 mg once daily for 5 days or until discharge if sooner) or usual standard of care alone (which includes dexamethasone 6 mg once daily for 10 days or until discharge if sooner). The primary outcome was 28-day mortality. On 11 May 2022, the independent Data Monitoring Committee recommended stopping recruitment of patients receiving no oxygen or simple oxygen only to this comparison due to safety concerns. We report the results for these participants only. Recruitment of patients receiving ventilatory support continues. The RECOVERY trial is registered with ISRCTN (50189673) and clinicaltrials.gov ( NCT04381936 ). Findings: Between 25 May 2021 and 12 May 2022, 1272 COVID-19 patients with hypoxia and receiving no oxygen (1%) or simple oxygen only (99%) were randomly allocated to receive usual care plus higher dose corticosteroids versus usual care alone (of whom 87% received low dose corticosteroids during the follow-up period). Of those randomised, 745 (59%) were in Asia, 512 (40%) in the UK and 15 (1%) in Africa. 248 (19%) had diabetes mellitus. Overall, 121 (18%) of 659 patients allocated to higher dose corticosteroids versus 75 (12%) of 613 patients allocated to usual care died within 28 days (rate ratio [RR] 1.56; 95% CI 1.18-2.06; p=0.0020). There was also an excess of pneumonia reported to be due to non-COVID infection (10% vs. 6%; absolute difference 3.7%; 95% CI 0.7-6.6) and an increase in hyperglycaemia requiring increased insulin dose (22% vs. 14%; absolute difference 7.4%; 95% CI 3.2-11.5). Interpretation: In patients hospitalised for COVID-19 with clinical hypoxia but requiring either no oxygen or simple oxygen only, higher dose corticosteroids significantly increased the risk of death compared to usual care, which included low dose corticosteroids. The RECOVERY trial continues to assess the effects of higher dose corticosteroids in patients hospitalised with COVID-19 who require non-invasive ventilation, invasive mechanical ventilation or extra-corporeal membrane oxygenation.

Apixaban following discharge in hospitalised adults with COVID-19: Preliminary results from a multicentre, open-label, randomised controlled platform clinical trial. (preprint)

Background The role of thromboprophylaxis in the post-acute phase of COVID-19 is uncertain due to conflicting results from randomised controlled trials and observational studies. We aimed to determine the effectiveness of post-hospital apixaban in reducing the rate of death and hospital readmission of hospitalised adults with COVID-19. Methods HEAL COVID is an adaptive randomised open label multicentre platform trial recruiting participants from National Health Service Hospitals in the United Kingdom. Here we report the preliminary results of apixaban comparison of HEAL-COVID. Participants with a hospital admission related to confirmed COVID-19 and an expected date of discharge in the subsequent five days were randomised to either apixaban 2.5 mg twice daily or standard care (no anticoagulation) for 14 days. The primary outcome was hospital free survival at 12 months obtained through routine data sources. The trial was prospectively registered with ISRCTN (15851697) and Clincialtrials.gov (NCT04801940). Findings Between 19 May 2021 and 21 November 2022, 402 participants from 109 sites were randomised to apixaban and 399 to standard care. Seven participants withdrew from the apixaban group and one from the standard care group. Analysis was undertaken on an intention-to-treat basis. The apixaban arm was stopped on the recommendation of the oversight committees following an interim analysis due to no indication of benefit. Of the 402 participants randomised to apixaban, 117 experienced death or rehospitalisation during a median follow-up of 344.5 days (IQR 125 to 365), and 123 participants receiving standard care experienced death or rehospitalisation during a median follow-up of 349 days (IQR 124 to 365). There was no statistical difference in the rate of death and rehospitalisation (HR: 0.96 99%CI 0.69-1.34; p=0.75). Three participants in the apixaban arm experienced clinically significant bleeding during treatment. Interpretation Fourteen days of post-hospital anticoagulation with the direct oral anticoagulant apixaban did not reduce the rate of death or rehospitalisation of adults hospitalised with COVID-19. These data do not support the use of prophylactic post-hospital anticoagulation in adults with COVID-19. Funding HEAL-COVID is funded by the National Institute for Health and Care Research [NIHR133788] and the NIHR Cambridge Biomedical Research Centre [ BRC-1215-20014*].

Variation inERAP2has opposing effects on severe respiratory infection and autoimmune disease (preprint)

Introduction ERAP2 is an aminopeptidase involved in immunological antigen presentation. Genotype data in human samples from before and after the Black Death, an epidemic due to Yersinia pestis , have marked changes in population allele frequency of the common single nucleotide polymorphism (SNP) rs2549794. This SNP in strong linkage disequilibrium with a key splicing SNP in ERAP2 (rs2248374) and this suggests that variation at ERAP2 may be relevant for protection from infection. rs2549794 is also associated with Crohn’s disease and findings imply balancing selection between infection and autoimmune disease at this locus. There have been no large-scale prospective case-control studies of variation at ERAP2 and infection. Methods This study aimed to explore the association between variation at ERAP2 and a) infection, b) autoimmune disease, and c) parental longevity as a proxy for lifespan. Genome Wide Association Studies (GWAS) of these outcomes were identified in contemporary cohorts (UK Biobank, FinnGen, and GenOMICC). Effect estimates were extracted for rs2549794 and rs2248374. Additionally, cis expression and protein quantitative trait loci (QTLs) for ERAP2 were used in Mendelian randomisation analyses. Results Across all cohorts, the T allele (minor allele frequency of 0.4-0.5) of rs2549794 showed evidence of association with respiratory infection (odds ratio; OR for pneumonia 1.03; 95% CI 1.01-1.05; p = 0.014). Effect estimates were larger in bacterial rather than viral infection and larger for more severe phenotypes (OR for critical care admission with pneumonia 1.08; 95% CI 1.02-1.14, p = 0.008, OR for death from pneumonia 1.07; 95% CI 1.01-1.12; p = 0.014). In contrast, opposing effects were identified for Crohn’s disease (OR 0.86; 95% CI 0.82-0.90, p = 8.6 × 10 −9 ) and type 1 diabetes (OR 0.95; 95% CI 0.90-0.99, p = 0.02). No strong evidence for association was identified for sepsis. Carriage of the T allele was associated with increased age of parental death (beta in Z-scored years across both parents age at death 0.01, 95% CI 0.004-0.017, p = 0.002). Similar results were identified for rs2248374. In Mendelian randomisation analyses, increasing transcription or protein levels of ERAP2 were strongly associated with protection from respiratory infection, with opposing effects identified on Crohn’s disease and type 1 diabetes. Increased expression of ERAP2 was associated with reduced parental longevity. Conclusions Variation at ERAP2 is associated with severe respiratory infection in modern societies, with an opposing association with Crohn’s disease and type 1 diabetes. These data support the hypothesis that changes in allele frequencies in ERAP2 observed at the time of the Black Death reflect protection from infection, and suggest ongoing balancing selection at this locus driven by autoimmune and infectious disease

Evaluation of S/F94 as a proxy for COVID-19 severity (preprint)

Optimising statistical power in early-stage trials and observational studies accelerates discovery and improves the reliability of results. Ideally, intermediate outcomes should be continuously distributed and lie on the causal pathway between an intervention and a definitive outcome such as mortality. In order to optimise power for an intermediate outcome in the RECOVERY trial, we devised and evaluated a modification to a simple, pragmatic measure of oxygenation function - the SaO2/FIO2 (S/F) ratio. We demonstrate that, because of the ceiling effect in oxyhaemoglobin saturation, S/F ceases to reflect pulmonary oxygenation function at high values of SaO2. Using synthetic and real data, we found that the correlation of S/F with a gold standard (PaO2/FIO2, P/F ratio) improved substantially when measurements with SaO2 > 0.94 are excluded(Spearman r, synthetic data: S/F: 0.31; S/F94: 0.85). We refer to this measure as S/F94. In order to test the underlying assumptions and validity of S/F94 as a predictor of a definitive outcome (mortality), we collected an observational dataset including over 39,000 hospitalised patients with COVID-19 in the ISARIC4C study. We first demonstrated that S/F94 is predictive of mortality in COVID-19. We then compared the sample sizes required for trials using different outcome measures (S/F94, the WHO ordinal scale, sustained improvement at day 28 and mortality at day 28) ensuring comparable effect sizes. The smallest sample size was needed when S/F94 on day 5 was used as an outcome measure. To facilitate future study design, we provide an online user interface to quantify realworld power for a range of outcomes and inclusion criteria, using a synthetic dataset retaining the population-level clinical associations in real data accrued in ISARIC4C https://isaric4c.net/endpoints. We demonstrated that S/F94 is superior to S/F as a measure of pulmonary oxygenation function and is an effective intermediate outcome measure in COVID-19. It is a simple and non-invasive measurement, representative of disease severity and provides greater statistical power to detect treatment differences than other intermediate endpoints.

Dimethyl fumarate in patients admitted to hospital with COVID-19 (RECOVERY): a randomised, controlled, open-label, platform trial (preprint)

Background: Dimethyl fumarate (DMF) is an anti-inflammatory drug that has been proposed as a treatment for patients hospitalised with COVID-19. Methods: This randomised, controlled, open-label platform trial (Randomised Evaluation of COVID-19 Therapy [RECOVERY]), is assessing multiple possible treatments in patients hospitalised for COVID-19. In this initial assessment of DMF, performed at 27 UK hospitals, eligible and consenting adults were randomly allocated (1:1) to either usual standard of care alone or usual standard of care plus DMF 120mg twice daily for 2 days followed by 240mg twice daily for 8 days, or until discharge if sooner. The primary outcome was clinical status on day 5 measured on a seven-point ordinal scale, assessed using a proportional odds model. Secondary outcomes were time to sustained improvement in clinical status, time to discharge, day 5 peripheral blood oxygenation, day 5 C-reactive protein, and improvement in day 10 clinical status. The trial is registered with ISRCTN (50189673) and clinicaltrials.gov (NCT04381936). Findings: Between 2 March 2021 and 18 November 2021, 713 patients were enrolled in the DMF evaluation, of whom 356 were randomly allocated to receive usual care plus DMF, and 357 to usual care alone. 95% of patients were receiving corticosteroids as part of routine care. There was no evidence of a beneficial effect of DMF on clinical status at day 5 (common odds ratio of unfavourable outcome 1.12; 95% CI 0.85-1.46; p=0.42). There was no significant effect of DMF on any secondary outcome. As expected, DMF caused flushing and gastrointestinal symptoms, each in around 6% of patients, but no new adverse effects were identified. Interpretation: In adults hospitalised with COVID-19, DMF was not associated with an improvement in clinical outcomes.

Nasal IgA wanes 9 months after hospitalisation with COVID-19 and is not induced by subsequent vaccination. (preprint)

Background Most studies of immunity to SARS-CoV-2 focus on circulating antibody, giving limited insights into mucosal defences that prevent viral replication and onward transmission. We studied nasal and plasma antibody responses one year after hospitalisation for COVID-19, including a period when SARS-CoV-2 vaccination was introduced. Methods Plasma and nasosorption samples were prospectively collected from 446 adults hospitalised for COVID-19 between February 2020 and March 2021 via the ISARIC4C and PHOSP-COVID consortia. IgA and IgG responses to NP and S of ancestral SARS-CoV-2, Delta and Omicron (BA.1) variants were measured by electrochemiluminescence and compared with plasma neutralisation data. Findings Strong and consistent nasal anti-NP and anti-S IgA responses were demonstrated, which remained elevated for nine months. Nasal and plasma anti-S IgG remained elevated for at least 12 months with high plasma neutralising titres against all variants. Of 180 with complete data, 160 were vaccinated between 6 and 12 months; coinciding with rises in nasal and plasma IgA and IgG anti-S titres for all SARS-CoV-2 variants, although the change in nasal IgA was minimal. Samples 12 months after admission showed no association between nasal IgA and plasma IgG responses, indicating that nasal IgA responses are distinct from those in plasma and minimally boosted by vaccination. Interpretation The decline in nasal IgA responses 9 months after infection and minimal impact of subsequent vaccination may explain the lack of long-lasting nasal defence against reinfection and the limited effects of vaccination on transmission. These findings highlight the need to develop vaccines that enhance nasal immunity.

Outcome of COVID-19 in hospitalised immunocompromised patients: an analysis of the WHO ISARIC CCP-UK prospective cohort study. (preprint)

Background Immunocompromised patients may be at higher risk of mortality if hospitalised with COVID-19 compared with immunocompetent patients. However, previous studies have been contradictory. We aimed to determine whether immunocompromised patients were at greater risk of in-hospital death, and how this risk changed over the pandemic. Methods We included patients >=19yrs with symptomatic community-acquired COVID-19 recruited to the ISARIC WHO Clinical Characterisation Protocol UK. We defined immunocompromise as: immunosuppressant medication preadmission, cancer treatment, organ transplant, HIV, or congenital immunodeficiency. We used logistic regression to compare the risk of death in both groups, adjusting for age, sex, deprivation, ethnicity, vaccination and co-morbidities. We used Bayesian logistic regression to explore mortality over time. Findings Between 17/01/2020 and 28/02/2022 we recruited 156,552 eligible patients, of whom 21,954 (14%) were immunocompromised. 29% (n=6,499) of immunocompromised and 21% (n=28,608) of immunocompetent patients died in hospital. The odds of inhospital mortality were elevated for immunocompromised patients (adjOR 1.44, 95% CI 1.39-1.50, p<0.001). As the pandemic progressed, in-hospital mortality reduced more slowly for immunocompromised patients than for immunocompetent patients. This was particularly evident with increasing age: the probability of the reduction in hospital mortality being less for immunocompromised patients aged 50-69yrs was 88% for men and 83% for women, and for those >80yrs was 99% for men, and 98% for women. Conclusions Immunocompromised patients remain at elevated risk of death from COVID-19. Targeted measures such as additional vaccine doses and monoclonal antibodies should be considered for this group.

Therapeutic potential of IL6R blockade for the treatment of sepsis and sepsis-related death: Findings from a Mendelian randomisation study (preprint)

Introduction Sepsis is characterised by dysregulated, life-threatening immune responses, which are thought to be driven by cytokines such as interleukin-6 (IL-6). Genetic variants in IL6R known to downregulate IL-6 signalling are associated with improved COVID-19 outcomes, a finding later confirmed in randomised trials of IL-6 receptor antagonists (IL6RA). We hypothesised that blockade of IL6R could also improve outcomes in sepsis. Methods We performed a Mendelian randomisation analysis using single nucleotide polymorphisms (SNPs) in and near IL6R to evaluate the likely causal effects of IL6R blockade on sepsis, sepsis severity, other infections, and COVID-19. We weighted SNPs by their effect on CRP and combined results across them in inverse variance weighted meta-analysis, proxying the effect of IL6RA. Our outcomes were measured in UK Biobank, FinnGen, the COVID-19 Host Genetics Initiative (HGI), and the GenOSept and GainS consortium. We performed several sensitivity analyses to test assumptions of our methods, including utilising variants around CRP in a similar analysis. Results In the UK Biobank cohort (N=485,825, including 11,643 with sepsis), IL6R blockade was associated with a decreased risk of sepsis (OR=0.80; 95% CI 0.66-0.96, per unit of natural log transformed CRP decrease). The size of this effect increased with severity, with larger effects on 28-day sepsis mortality (OR=0.74; 95% CI 0.38-0.70); critical care admission with sepsis (OR=0.48, 95% CI 0.30-0.78) and critical care death with sepsis (OR=0.37, 95% CI 0.14 - 0.98) Similar associations were seen with severe respiratory infection: OR for pneumonia in critical care 0.69 (95% CI 0.49 - 0.97) and for sepsis survival in critical care (OR=0.22; 95% CI 0.04- 1.31) in the GainS and GenOSept consortium. We also confirm the previously reported protective effect of IL6R blockade on severe COVID-19 (OR=0.69, 95% 0.57 - 0.84) in the COVID-19 HGI, which was of similar magnitude to that seen in sepsis. Sensitivity analyses did not alter our primary results. Conclusions IL6R blockade is causally associated with reduced incidence of sepsis, sepsis related critical care admission, and sepsis related mortality. These effects are comparable in size to the effect seen in severe COVID-19, where IL-6 receptor antagonists were shown to improve survival. This data suggests a randomised trial of IL-6 receptor antagonists in sepsis should be considered.

An international observational study to assess the impact of the Omicron variant emergence on the clinical epidemiology of COVID-19 in hospitalised patients (preprint)

BackgroundWhilst timely clinical characterisation of infections caused by novel SARS-CoV-2 variants is necessary for evidence-based policy response, individual-level data on infecting variants are typically only available for a minority of patients and settings. MethodsHere, we propose an innovative approach to study changes in COVID-19 hospital presentation and outcomes after the Omicron variant emergence using publicly available population-level data on variant relative frequency to infer SARS-CoV-2 variants likely responsible for clinical cases. We apply this method to data collected by a large international clinical consortium before and after the emergence of the Omicron variant in different countries. ResultsOur analysis, that includes more than 100,000 patients from 28 countries, suggests that in many settings patients hospitalised with Omicron variant infection less often presented with commonly reported symptoms compared to patients infected with pre-Omicron variants. Patients with COVID-19 admitted to hospital after Omicron variant emergence had lower mortality compared to patients admitted during the period when Omicron variant was responsible for only a minority of infections (odds ratio in a mixed-effects logistic regression adjusted for likely confounders, 0.67 [95% confidence interval 0.61 - 0.75]). Qualitatively similar findings were observed in sensitivity analyses with different assumptions on population-level Omicron variant relative frequencies, and in analyses using available individual-level data on infecting variant for a subset of the study population. ConclusionsAlthough clinical studies with matching viral genomic information should remain a priority, our approach combining publicly available data on variant frequency and a multi-country clinical characterisation dataset with more than 100,000 records allowed analysis of data from a wide range of settings and novel insights on real-world heterogeneity of COVID-19 presentation and clinical outcome.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative (preprint)

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,048 severe disease cases and 571,009 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p=5.41x10-7). These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

GWAS and meta-analysis identifies multiple new genetic mechanisms underlying severe Covid-19. (preprint)

Pulmonary inflammation drives critical illness in Covid-19, creating a clinically homogeneous extreme phenotype, which we have previously shown to be highly efficient for discovery of genetic associations. Despite the advanced stage of illness, we have found that immunomodulatory therapies have strong beneficial effects in this group. Further genetic discoveries may identify additional therapeutic targets to modulate severe disease. In this new data release from the GenOMICC (Genetics Of Mortality in Critical Care) study we include new microarray genotyping data from additional critically-ill cases in the UK and Brazil, together with cohorts of severe Covid-19 from the ISARIC4C and SCOURGE studies, and meta-analysis with previously-reported data. We find an additional \numconvincingnew new genetic associations. Many are in potentially druggable targets, in inflammatory signalling (JAK1, PDE4A), monocyte-macrophage differentiation (CSF2), immunometabolism (SLC2A5, AK5), and host factors required for viral entry and replication (TMPRSS2, RAB2A). As with our previous work, these results provide tractable therapeutic targets for modulation of harmful host-mediated inflammation in Covid-19.

Baricitinib in patients admitted to hospital with COVID-19 (RECOVERY): a randomised, controlled, open-label, platform trial and updated meta-analysis (preprint)

Background: We evaluated the use of baricitinib, a Janus kinase (JAK) 1/2 inhibitor, for the treatment of patients admitted to hospital because of COVID-19. Methods: This randomised, controlled, open-label platform trial (Randomised Evaluation of COVID-19 Therapy [RECOVERY]), is assessing multiple possible treatments in patients hospitalised for COVID-19. Eligible and consenting patients were randomly allocated (1:1) to either usual standard of care alone (usual care group) or usual care plus baricitinib 4 mg once daily by mouth for 10 days or until discharge if sooner (baricitinib group). The primary outcome was 28-day mortality assessed in the intention-to-treat population. A meta-analysis was conducted that included the results from the RECOVERY trial and all previous randomised controlled trials of baricitinib or other JAK inhibitor in patients hospitalised with COVID-19. Findings: Between 2 February 2021 and 29 December 2021, 8156 patients were randomly allocated to receive usual care plus baricitinib versus usual care alone. At randomisation, 95% of patients were receiving corticosteroids and 23% receiving tocilizumab (with planned use within the next 24 hours recorded for a further 9%). Overall, 513 (12%) of 4148 patients allocated to baricitinib versus 546 (14%) of 4008 patients allocated to usual care died within 28 days (age-adjusted rate ratio 0.87; 95% CI 0.77-0.98; p=0.026). This 13% proportional reduction in mortality was somewhat smaller than that seen in a meta-analysis of 8 previous trials of a JAK inhibitor (involving 3732 patients and 425 deaths) in which allocation to a JAK inhibitor was associated with a 43% proportional reduction in mortality (rate ratio 0.57; 95% CI 0.45-0.72). Including the results from RECOVERY into an updated meta-analysis of all 9 completed trials (involving 11,888 randomised patients and 1484 deaths) allocation to baricitinib or other JAK inhibitor was associated with a 20% proportional reduction in mortality (rate ratio 0.80; 95% CI 0.71-0.89; p<0.001). In RECOVERY, there was no significant excess in death or infection due to non-COVID-19 causes and no excess of thrombosis, or other safety outcomes. Interpretation: In patients hospitalised for COVID-19, baricitinib significantly reduced the risk of death but the size of benefit was somewhat smaller than that suggested by previous trials. The total randomised evidence to date suggests that JAK inhibitors (chiefly baricitinib) reduce mortality in patients hospitalised for COVID-19 by about one-fifth.

Systematic comparison of ranking aggregation methods for gene lists in experimental results (preprint)

A common experimental output in biomedical science is a list of genes implicated in a given biological process or disease. The results of a group of studies answering the same, or similar, questions can be combined by meta-analysis to find a consensus or a more reliable answer. Ranking aggregation methods can be used to combine gene lists from various sources in meta-analyses. Evaluating a ranking aggregation method on a specific type of dataset before using it is required to support the reliability of the result since the property of a dataset can influence the performance of an algorithm. Evaluation of aggregation methods is usually based on a simulated database especially for the algorithms designed for gene lists because of the lack of a known truth for real data. However, simulated datasets tend to be too small compared to experimental data and neglect key features, including heterogeneity of quality, relevance and the inclusion of unranked lists. In this study, a group of existing methods and their variations which are suitable for meta-analysis of gene lists are compared using simulated and real data. Simulated data was used to explore the performance of the aggregation methods as a function of emulating the common scenarios of real genomics data, with various heterogeneity of quality, noise level, and a mix of unranked and ranked data using 20000 possible entities. In addition to the evaluation with simulated data, a comparison using real genomic data on the SARS-CoV-2 virus, cancer (NSCLC), and bacteria (macrophage apoptosis) was performed. We summarise our evaluation results in terms of a simple flowchart to select a ranking aggregation method for genomics data.

Rapid selection of P323L in the SARS-CoV-2 polymerase (NSP12) in humans and non-human primate models and confers a large plaque phenotype (preprint)

The mutational landscape of SARS-CoV-2 varies at both the dominant viral genome sequence and minor genomic variant population. An early change associated with transmissibility was the D614G substitution in the spike protein. This appeared to be accompanied by a P323L substitution in the viral polymerase (NSP12), but this latter change was not under strong selective pressure. Investigation of P323L/D614G changes in the human population showed rapid emergence during the containment phase and early surge phase of wave 1 in the UK. This rapid substitution was from minor genomic variants to become part of the dominant viral genome sequence. A rapid emergence of 323L but not 614G was observed in a non-human primate model of COVID-19 using a starting virus with P323 and D614 in the dominant genome sequence and 323L and 614G in the minor variant population. In cell culture, a recombinant virus with 323L in NSP12 had a larger plaque size than the same recombinant virus with P323. These data suggest that it may be possible to predict the emergence of a new variant based on tracking the distribution and frequency of minor variant genomes at a population level, rather than just focusing on providing information on the dominant viral genome sequence e.g., consensus level reporting. The ability to predict an emerging variant of SARS-CoV-2 in the global landscape may aid in the evaluation of medical countermeasures and non-pharmaceutical interventions.

Comparison of children and young people admitted with SARS-CoV-2 across the UK in the first and second pandemic waves: prospective multicentre observational cohort study (preprint)

BackgroundChildren and young people (CYP) were less affected than adults in the first wave of SARS-CoV-2 in the UK. We test the hypothesis that clinical characteristics of hospitalized CYP with SARS-CoV-2 in the UK second wave would differ from the first due to the combined impact of the alpha variant, school reopening and relaxation of shielding. MethodsPatients <19 years hospitalised in the UK with clinician-reported SARS-CoV-2 were enrolled in a prospective multicentre observational cohort study between 17th January 2020 and 31st January 2021. Minimum follow up time was two weeks. Clinical characteristics were compared between the first (W1) and second wave (W2) of infections. Findings2044 CYP aged <19 years were reported from 187 hospitals. 427/2044 (20.6%) had asymptomatic/incidental SARS-CoV-2 infection and were excluded from main analysis. 16.0% (248/1548) of symptomatic CYP were admitted to critical care and 0.8% (12/1504) died. 5.6% (91/1617) of symptomatic CYP had Multisystem Inflammatory Syndrome in Children (MIS-C). Patients in W2 were significantly older (median age 6.5 years, IQR 0.3-14.9) than W1 (4.0 (0.4-13.6, p 0.015). Fever was more common in W1, otherwise presenting symptoms and comorbidities were similar across waves. After excluding CYP with MIS-C, patients in W2 had lower PEWS at presentation, lower antibiotic use and less respiratory and cardiovascular support compared to W1. There was no change in the proportion of CYP admitted to critical care between W1 and W2. 58.0% (938/1617) of symptomatic CYP had no reported comorbidity. Patients without co-morbidities were younger (42.4%, 398/938, <1 year old), had lower Paediatric Early Warning Scores (PEWS) at presentation, shorter length of hospital stay and received less respiratory support. MIS-C was responsible for a large proportion of critical care admissions, invasive and non-invasive ventilatory support, inotrope and intravenous corticosteroid use in CYP without comorbidities. InterpretationSevere disease in CYP admitted with symptomatic SARS-CoV-2 in the UK remains rare. One in five CYP in this cohort had asymptomatic/incidental SARS-CoV-2 infection. We found no evidence of increased disease severity in W2 compared with W1. FundingShort form: National Institute for Health Research, UK Medical Research Council, Wellcome Trust, Department for International Development and the Bill and Melinda Gates Foundation. Long form: This work is supported by grants from the National Institute for Health Research (award CO-CIN-01) and the Medical Research Council (grant MC_PC_19059) and by the National Institute for Health Research Health Protection Research Unit (NIHR HPRU) in Emerging and Zoonotic Infections at University of Liverpool in partnership with Public Health England (PHE), in collaboration with Liverpool School of Tropical Medicine and the University of Oxford (NIHR award 200907), Wellcome Trust and Department for International Development (215091/Z/18/Z), and the Bill and Melinda Gates Foundation (OPP1209135). Liverpool Experimental Cancer Medicine Centre provided infrastructure support for this research (grant reference: C18616/A25153). JSN-V-T is seconded to the Department of Health and Social Care, England (DHSC). The views expressed are those of the authors and not necessarily those of the DHSC, DID, NIHR, MRC, Wellcome Trust, or PHE.

Common, low-frequency, rare, and ultra-rare coding variants contribute to COVID-19 severity. (preprint)

The combined impact of common and rare exonic variants in COVID-19 host genetics is currently insufficiently understood. Here, common and rare variants from whole exome sequencing data of about 4,000 SARS-CoV-2-positive individuals were used to define an interpretable machine learning model for predicting COVID-19 severity. Firstly, variants were converted into separate sets of Boolean features, depending on the absence or the presence of variants in each gene. An ensemble of LASSO logistic regression models was used to identify the most informative Boolean features with respect to the genetic bases of severity. The Boolean features selected by these logistic models were combined into an Integrated PolyGenic Score that offers a synthetic and interpretable index for describing the contribution of host genetics in COVID-19 severity, as demonstrated through testing in several independent cohorts. Selected features belong to ultra-rare, rare, low-frequency, and common variants, including those in linkage disequilibrium with known GWAS loci. Noteworthly, around one quarter of the selected genes are sex-specific. Pathway analysis of the selected genes associated with COVID-19 severity reflected the multi-organ nature of the disease. The proposed model might provide useful information for developing diagnostics and therapeutics, while also being able to guide bedside disease management.

Whole genome sequencing identifies multiple loci for critical illness caused by COVID-19 (preprint)

Critical illness in COVID-19 is caused by inflammatory lung injury, mediated by the host immune system. We and others have shown that host genetic variation influences the development of illness requiring critical care1 or hospitalisation2;3;4 following SARS-Co-V2 infection. The GenOMICC (Genetics of Mortality in Critical Care) study is designed to compare genetic variants in critically-ill cases with population controls in order to find underlying disease mechanisms. Here, we use whole genome sequencing and statistical fine mapping in 7,491 critically-ill cases compared with 48,400 population controls to discover and replicate 22 independent variants that significantly predispose to life-threatening COVID-19. We identify 15 new independent associations with critical COVID-19, including variants within genes involved in interferon signalling (IL10RB, PLSCR1), leucocyte differentiation (BCL11A), and blood type antigen secretor status (FUT2). Using transcriptome-wide association and colocalisation to infer the effect of gene expression on disease severity, we find evidence implicating expression of multiple genes, including reduced expression of a membrane flippase (ATP11A), and increased mucin expression (MUC1), in critical disease. We show that comparison between critically-ill cases and population controls is highly efficient for genetic association analysis and enables detection of therapeutically-relevant mechanisms of disease. Therapeutic predictions arising from these findings require testing in clinical trials.

Pandemic, epidemic, endemic: B cell repertoire analysis reveals unique anti-viral responses to SARS-CoV-2, Ebola and Respiratory Syncytial Virus (preprint)

Immunoglobulin gene heterogeneity reflects the diversity and focus of the humoral immune response towards different infections, enabling inference of B cell development processes. Detailed compositional and lineage analysis of long read IGH repertoire sequencing, combining examples of pandemic, epidemic and endemic viral infections with control and vaccination samples, demonstrates general responses including increased use of IGHV4-39 in both EBOV and COVID-19 infection cohorts. We also show unique characteristics absent in RSV infection or yellow fever vaccine samples: EBOV survivors show unprecedented high levels of class switching events while COVID-19 repertoires from acute disease appear underdeveloped. Despite the high levels of clonal expansion in COVID-19 IgG1 repertoires there is a striking lack of evidence of germinal centre mutation and selection. Given the differences in COVID-19 morbidity and mortality with age, it is also pertinent that we find significant differences in repertoire characteristics between young and old patients. Our data supports the hypothesis that a primary viral challenge can result in a strong but immature humoral response where failures in selection of the repertoire risks off-target effects.